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Recent developements in QCD at finite temperature are reviewed. Particular emphasis is laid on results 
stemming from simulations which involve quarks. 



1. INTRODUCTION 

Lattice investigations of QCD at finite temper- 
ature have contributed considerably to the cur- 
rent understanding of the transition from the 
hadronic state of matter to the quark gluon 
plasma and of the physics of the plasma phase. 
Mainly due to the ever-present limitations of com- 
putational power many analyses have been car- 
ried out in the pure gauge sector of QCD such 
that bulk properties of gluons at finite tempera- 
ture can be regarded as solved: the system has 
a well-established first-order transition ||^, the 
equation of state is known in the continuum limit 
|2| and the critical temperature in the contin- 
uum limit has been determined with only a few 
percent uncertainty [HJ|]. Clearly, more detailed 
questions like e.g. the nature of excitations in the 
plasma deserve further work, also in the quenched 
approximation. Yet, the emphasis of recent re- 
search has shifted towards studies of full QCD in- 
cluding staggered as well as Wilson quarks. These 
studies consistently have lead to an estimate of 
the critical temperature of order 150 MeV for 2 
flavors so far. This value is considerably lower 
than the quenched number of = 270(5)MeV. 
Through the relation T — l/{aNr) of the temper- 
ature T to lattice spacing a and temporal extent 
Nr of the lattice, dynamical fermion simulations 
in the vicinity of the transition are, at a given Nt, 
carried out at considerably larger lattice spac- 
ings. Already for that reason, with the standard 
discretizations, extrapolations to the contiunuum 
limit will be more difficult than in the quenched 
case. Therefore, the search for improved actions 
has received much attention recently in the con- 
text of finite temperature QCD. 



This review attempts to summarize the devel- 
opements in finite temperature lattice QCD since 
last year's conference as reviewed in I will 
concentrate on problems and results in simula- 
tions with full QCD. Section 2 summarizes esti- 
mates of the critical temperature. In section 3 
the present status of knowledge about the nature 
of the chiral transition is discussed, mainly for 
two flavors staggered as well as Wilson fermions. 
Section 4 describes studies of energy density and 
pressure in the high temperature limit while sec- 
tion 5 reviews a few recent results on screening 
lengths and masses. This year's results at finite 
density are summarized in section 6, conclusions 
are given in section 7. 

2. CRITICAL TEMPERATURE 

One of the basic quantities to be derived from 
finite temperature lattice QCD is the value of 
the critical temperature. In order to prepare 
the stage for the subsequent discussion of re- 
sults from dynamical fermion simulations, Fig- 
ure summarizes the current status of analy- 
sis in the quenched approximation Figure |^ 
shows the ratio Tc/^/a where a is the string 
tension extracted from the static quark poten- 
tial for various actions. For all data points the 
value of the critical coupling has been extrapo- 
lated to its infinite (spatial) volume limit at which 
then the string tension was determined. The 
lowest set of data points originates from simula- 
tions with the standard Wilson gauge action . 
A quadratic extrapolation in the lattice spacing 
to the continuum limit gives Tc/^/a = 0.631(2). 
Symanzik-improved actions show a much weaker 
cut-off dependence and are in agreement with the 
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continuum extrapolation of the standard action 
Likewise, the resuhs 1^ from Iwasaki's RG- 
improved action are consistent with a constant 
behavior in a, but they dehver a value of the crit- 
ical temperature of Tcj \/a — 0.656(4) which is 
about 3% higher than the number from the stan- 
dard action. Since the procedure to extract the 
string tension has not been the same for the two 
numbers, one might suspect that the difference in 
the quoted value for Tc, = 276(2)MeV versus 
266(2) MeV, is mainly due to differences in the 
analysis of the static quark potential Q rather 
than to differences in the improvement scheme. 
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Figure 1. The quenched critical temperature in 
units of the square root of the string tension for 
various gauge actions versus the the lattice spac- 
ing squared. 



The current situation with dynamical fermions 
is depicted in Figure I The plot summarizes data 
from simulations with 2 flavors of quarks, stag- 
gered fermions at A^^- = 4 and 6 (H,^ as well as 
improved Wilson fermions at Nr = 4 plus 
Np = 4 staggered results obtained from Nr ~ 4 
lattices with an improved action in addition 
to an old number from N^- — 8 and a standard 
action. Compared to the equivalent quenched 
plot. Figure shows that the lattice spacings at 



which Tc has been determined so far are consider- 
ably larger than in pure gauge theory simulations. 
Moreover, the investigations have not been car- 
ried out at the physical quark masses. The arrow 
in Figure ^ indicates that at fixed the tran- 
sition takes place at larger lattice spacings when 
the quark mass is decreased. Thus the critical 
temperature is decreasing when the quark mass 
is lowered. 
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Figure 2. The critical temperature in units of the 
square root of the string tension for dynamical 
fermions versus the square of the lattice spacing 
(for further explanations see text). 



The same data is shown again as a function of 
the pseudoscalar Goldstone boson to vector me- 
son mass ratio {Mps /MvY in Figure ||. (The 
meson masses at the appropriate critical values 
of the coupling in parts have been obtained by 
means of a phenomenological interpolation for- 
mula |l|].) Here now, at fixed Nr, smaller lattice 
spacings are to the right of the figure. The Nr = 4 
staggered data indicates that over ^/a tends 
to lower values as the quark mass is decreased. 
The same trend is observed for the Wilson im- 
proved results, although at larger [Adps/MvY 
ratio. On the other hand, the Nr — 6 data point 
seems to indicate that, at a given quark mass, de- 
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Figure 3. The critical temperature in units of the 
square root of the string tension for dynamical 
fermions plotted versus [Mps/MvY ■ The point 
at (Mps /MvY = 1 is the 7V^ = 4 quenched value. 
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Figure 4. The critical temperature in units of the 
vector meson mass for staggered fermions plotted 
versus {Mps/Mvf- 



creasing the lattice spacing increases Tc/ ^/a only 
slightly. At the moment, one would therefore esti- 
mate a physical value for the critical temperature 
of Tc/v^<0.4 or Tc<170MeV. 

The critical temperature has also been esti- 
mated from the ratio to the vector meson mass. 
In this case one ought to go (close) to the chiral 
limit in order to extract a physical number. In 
the case of using the string tension to set the scale 
one might argue that the string tension is consid- 
erably less affected by the quark mass. Figure ^ 
shows Tel My for Np = 2 staggered fermions 
gjl^^jll, plotted as function of {Mps/Myf. 
As the quark mass is decreased this ratio rises. 
Recall that the infinite quark mass, quenched 
data point corresponds to Mps/My = 1 and 
Tc/Mv = 0. As the lattice spacing is decreased, 
Tc/My stays remarkably constant. Extrapolat- 
ing the A'^^ = 4 data to the chiral limit suggests 
a value of Tc/My ^ 0.2 or ~ ISOMeV. Note 
that this value disagrees somewhat with the num- 
ber extracted from the string tension. 

The corresponding data for dynamical Wil- 
son quarks |p|,p7|-p^ are given in Figure ||. De- 
spite the known problems with standard Wilson 



ferniion thermodynamics, at least the results for 
Nr — 6 and 8 as well as the first data with im- 
proved Wilson fermions Q are not in disagree- 
ment with the staggered data. 

3. PHASE TRANSITION 

3.1. Staggered Np — 2: critical behavior 

The theoretical expectations on the scaling be- 
havior of the theory at the chiral transition are 
based on the cr model in three dimensions. For the 
case of two light flavors, if the transition is second 
order, it is expected to show scaling behavior with 
SU{2) X SU{2) ~ 0(4) exponents. On the other 
hand, if the anomalous symmetry were ef- 

fectively restored, the relevant symmetry group 
would be UAil) x 5';7(2) x S'C/(2) ~ C'(2£x 0(4) 
and the transition could be first order ||2^ . 

It has been attempted to analyze the critical 
behavior of 2 flavor staggered QCD by studying 
the scaling behavior of various quantities and de- 
termining critical exponents |21|. These scaling 
relations are derived from the scaling of the sin- 
gular part of the free energy density, 

fit, h) = -^\nz = h-^fibyn, by^h) (i) 
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Figure 5. The critical temperature in units of 
the vector meson mass for Wilson fermions versus 

{Mps/Mvf. 



Here T and V denote temperature and volume, 
Z is the partition function, t is the reduced tem- 
perature, t — (T — Tc)/Tc, and h is the symmetry 
breaking field, h — m/T where m is the quark 
mass. The scaling factor 6 is arbitrary and can be 
chosen appropriately. In the vicinity of the crit- 
ical point thermodynamic quantities should be 
governed by the thermal {yt) and the magnetic 
{yh) critical exponent. In the staggered version 
of lattice regularized QCD, for the dimensionless 
couplings t and h one uses 

6 6 
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h = maN^ 



(2) 



where 5c (0) denotes the critical coupling on a lat- 
tice with fixed temporal extent in the limit of van- 
ishing quark mass. At non- vanishing quark mass, 
a pseudo-critical coupling gdm) is defined as the 
location of a peak in e.g. the Polyakov loop sus- 
ceptibility. 

The status of investigations of the chiral tran- 
sition with two light staggered quarks has been 
commented upon at last year's conference |^]. 
Since then, JLQCD as well as the Bielefeld 
group |23] have about finalized their analyses. A 



new analysis was presented by C. DeTar at this 
conference 24 1 . 

The quantities analyzed by JLQCD and the 
Bielefeld group are various susceptibilities, in par- 
ticular the magnetic or chiral susceptibility 



V ^ dm 

i—l 



and the thermal susceptibility 



T 



Xt 



V ^ dm.,d{l/T) 



(3) 



(4) 



Assuming that the free energy is dominated by 
its singular part, Eq. ([l]) then leads to the scaling 
predictions for the peak heights of the suscepti- 
bilities at the line of pseudo-critical couplings 



peak 
Am 
peak 

Xt 



(5) 



where the exponents are given hy z„i = 2 — 1 jy^ 
and zt = (yt — + 1. The pseudo-critical line 
itself is expected to follow 
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gl(jn) g2(o) 



cm 



(6) 



with Zg = yt/yh- The values of these exponents 
for various symmetries are given in table |l|. 
At finite lattice spacing the exact chiral symmetry 
of the staggered fermion action is C/(l) — 0(2). 
However, sufficiently close to the continuum limit 
one expects 0(4) exponents. The possibility of 
mean-field (MF) exponents arbitrarily close to 
the transition has been raised by . 

Earlier investigations of the exponents on small 
lattices (8'^ x 4) had observed partial agreement 
with 0(4) scaling ||2l[. These studies have been 
repeated on larger spatial volumes, L — 12, 16. 
In addition to the quark mass values 0.02, 0.0375 
and 0.075 JLQCD also ran at m = 0.01. The vol- 
ume dependence of the chiral peak susceptibility 
is shown in Figure |^, similar results are available 
for the other quantities. For m > 0.02, the sus- 
ceptibility rises when the volume is increased from 
8'^ to 12'^, but then stays approximately constant. 
Thus, a phase transition does not occur in this 
mass range, in agreement with earlier claims p7|] . 



5 



Table 1 

Critical exponents for 0(2),0(4) and mean field 
(MF). The 2 flavor QCD results are given sep- 
arately for each spatial lattice size, with upper 
values denoting the JLQCD and the lower ones 
the Bielefeld group numbers. 

0(2) 0(4) MF L^8 L=12 L=16 
Vg 0.60 0.54 2/3 0.70(11) 0.74(6) 0.64(5) 
0.63(6) 

0.79 0.79 2/3 0.70(4) 0.99(8) 1.03(9) 
0.84(5) 1.06(7) 0.93(8) 
zt 0.39 0.34 1/3 0.47(5) 0.81(9) 0.83(12) 
0.63(7) 0.94(12) 0.85(12) 



Table y. For Zg, within two standard deviations 
agreement with all three predictions is obtained. 
For the other two exponents, both groups consis- 
tently observe a drastic change when the volume 
is increased from L — 8 to L — 12, 16. While 
for the small volume the value for z^ is in rough 
agreement with 0{2) and 0(4), the results from 
L = 12 and 16 do not agree with any of the 
predicted numbers. Indeed, the observed value 
z„i ~ 1 would be expected for a first order transi- 
tion. The thermal exponent zt is larger than any 
of the predictions for all volumes. 



At m — 0.01 the linear increase in the peak height 
as the volume is enlarged continues up to i = 16. 
As such, this observation could suggest a first or- 
der transition. JLQCD however have studied the 
volume dependence of a double-peak structure in 
the distribution of the chiral order parameter and 
conclude that a first order transition is likely to 
be absent I2I. 
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Figure 6. Volume dependence of the chiral sus- 
ceptibility at peak. 



The quark mass dependence of Xm'^^ is shown 
in Figure 0, together with fits to the expected 
scaling behavior, Eq. (^). The resulting values 
for the critical exponents are also summarized in 



Figure 7. Mass dependence of the chiral suscep- 
tibility at peak. The upper 4 lines are fit results 
with Eq. (^ to the L — 12 and 16 data while the 
lowest line shows the slope of 0(4) scaling. 



Another way to study the scaling behavior is 
to compute the (magnetic) equation of state |p8|| 



(7) 



and compare it with the scaling function 4> as de- 
termined from a parametrization of 0(4) simula- 
tion results [p9| , after adjusting two non- universal 
normalization constants. This has been done by 
the MILC Collaboration H for A^^ = 4, 6, 8 and 
12. The result for Nr = 4 is shown in Fig- 
ure 1^. While the data for the larger quark masses 
and smaller volumes are compatible with 0(4), 
the new data at smaller quark masses and larger 
lattice extent again show drastic disagreement. 
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Figure 8. The magnetic equation of state, Eq. (|^), 
at A^r = 4 The data at quark masses 0.0375, 
0.02 and 0.025 comes from L = 8 lattices while at 
0.0125 and 0.008 lattice extents up to L = 24 have 
been used. The line is the 0(4) scaling prediction. 
It can be moved horizontally as well as vertically 
by adjusting two free normalization constants. 



(MILC also has looked at one of the exponents, 
with the result zt = 1-3(2).) When is in- 
creased, thus going to smaller lattice spacings, the 
agreement becomes increasingly better p^ , but 
it should be remarked that the data at Nr — 12 
originates from quark mass values m/T ~ 0.1 
which are of about the same size as the larger 
quark masses used at A^t- =4. Also, even at large 
spacing one would expect 0(2) behavior which 
is indistinguishable from 0(4) with the current 
precision of the data. 

At the moment there is no convincing explana- 
tion for these discrepancies at hand. In view of 
the results presented in section 3.3 obtained with 
an improved gauge action and Wilson fermions 
one might speculate that at strong coupling and 
for the standard action the relation Eq. ^ be- 
tween the QCD parameters and the thermody- 
namic variables as they enter the singular part of 
the free energy is strongly distorted. More stud- 
ies at weaker coupling or with improved actions 



would be needed to solve this important question. 

3.2. Ua{1) restoration 

The nature of the chiral transition for two fla- 
vors is strongly affected by the realization of the 
[/yi(l) symmetry ||2^. At very high temperatures 
topologically non-trivial configurations are sup- 
pressed, thus leading to the effective restoration 
of the symmetry despite the anomaly. For 2 light 
quark flavors the effective restoration of Ua{1) is 
reflected in the degeneracy of the pion and the 
isovector-scalar ao (S) mass Q . This degeneracy 
can also be detected by comparing the correlation 
functions 

u = j d^x{{^{x)n{Q)) ~ (ao(a;)ao(0))) (8) 

If Ua{1) is restored this quantity should vanish in 
the chiral limit. The (generalized) susceptibility 
has been looked at by various groups pl| , ^ , p3| . 
At finite quark mass uo is dropping across the chi- 
ral transition, but stays non-zero above Tc- Fig- 
ure H shows the latest Columbia data for this 
quantity, together with various extrapolations to 
the chiral limit. 

In the continuum, the susceptibility lo is ex- 
pected to be an analytic and, for Np = 2, even 
function in the quark mass. Indeed, fits with a 
quadratic m dependence work and lead to a fi- 
nite intercept in the chiral limit. However, the 
data look strikingly linear and fitting them with 
a linear ansatz results in a vanishing of the sus- 
ceptibility at m = 0. At finite lattice spacing, due 
to zero-mode shifts and taking the square root of 
the determinant the approach towards the chiral 
limit is not so clear |Q. Therefore one should 
continue to study the quark mass dependence at 
even smaller quark masses as well as at smaller 
lattice spacings. 

The approach chosen in is to determine 

screening masses. Above the critical temperature, 
the difference between tt and ag mass drops con- 
siderably, but a non-degeneracy remains at finite 
quark mass, thus confirming the findings origi- 
nating from the analysis of the susceptibility. In 
order to address the problem of the chiral limit 
from a different angle, ref. ||3j] also computed the 
lowest eigenvalues A and corresponding eigenvec- 
tors f/'A of the fermion matrix. In the continuum, 
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Figure 9. The quantity Eq. (g|), which mea- 
sures the breaking of the \J a(X) symmetry plot- 
ted versus the quark mass (lowest data) together 
with various fits js^. Also shown are the results 
for the integrated pion correlator, xp, and the 
chiral condensate. The data were obtained on a 
16^ X 4 lattice at fixed /? slightly above /3c- 



in the phase symmetric with respect to the ax- 
ial SU{2), the chiral limit of iv is given by the 
zero-modes, 



(E 

A=0 



iX + m ' 



(9) 



In |34| it is then verified that ui obtained in the 
standard fashion is saturated by the contribution 
from low eigenmodes thus supporting that contin- 
uum arguments seem to hold at lattices spacings 
of ©(O.lSfm) {Nr = 8). The contribution of zero 
modes in Eq. (|^) is to leading order quark mass 
independent. If this could be verified e.g. by run- 
ning at more quark mass values the results of 
would indicate that the Ua{^) symmetry is not 
restored at the chiral transition. 



3.3. Wilson fermions Np = 2 

The phase diagram of QCD at finite temper- 
ature with 2 flavors of Wilson quarks has been 
clarified in ||35|| and has been explained in great 
detail in last year's review At finite Nr, the 
line Kc(/3) defined through the vanishing of the 
pion mass starts off at 1/4 at /3 = and extends 
to Kc ~ 0.22 at about /3 ~ 4.0 where it bends 
backwards again to the region of stronger cou- 
phngs (see also Q). On the other hand, coming 
from the confined phase, at the thermal line Kt{j3) 
where the Polyakov loop develops a non- vanishing 
expectation value the pion mass increases rapidly 
due to the approximate restoration of chiral sym- 
metry. Only in the region where the thermal line 
is close to Kc does the theory have a pion with 
a small mass. Thus, the chiral transition can 
only be explored in that region. Unfortunately, 
this region is at strong coupling for iVi- = 4 and 
moves towards smaller coupling only very slowly 
with increasing temporal extent of the lattice , 
rendering a study of the transition in the vicin- 
ity of continuum physics prohibitively expensive. 
For that reason and for the well-known patholo- 
gies several groups have started to work with 
improved actions. 

Iwasaki et al. simulated with the standard Wil- 
son fermion action but on RG-improved glue [^8| . 
Qualitatively, the phase diagram is very similar 
to the standard one so that small pion masses 
again are obtained in the vicinity of the finite tem- 
perature Kc cusp. In addition to the phase dia- 
gram the group has also investigated the magnetic 
equation of state, Eq. (Q). For Wilson fermions 
quark mass and chiral order parameter have to be 
obtained from chiral Ward identities This 
involves renormalization constants for which the 
lowest order perturbative values have been used 
in 1^^. The results, including new data at low /3 
are shown in Figure |l^. The agreement with the 
0(4) scaling curve is remarkable. The analysis 
was carried out on lattices of size 8^ x 4 and with 
mainly not very small quark masses. It would 
be very interesting to continue the investigation 
on larger lattices and with more data at smaller 
quark masses. 

The MILC Collaboration has analyzed the fi- 
nite temperature transition with the Symanzik- 
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scaling analysis with 0(4) exponents 

0(4) beta_ct=1 .356 chi"2/df=0.71 




t/h»1/(beta*clelta) 

Figure 10. The magnetic equation of state, 
Eq. (0), with two flavors of standard Wilson 
fermions on improved glue [p8[ . 



Liischer-Weisz gluon action in combination with 
the Sheikholeslami-Wohlert fermion action, both 
tadpole improved They obtain a much 

smoother cross-over behavior as compared with 
the standard actions which suggests that the 
pathologies of the latter are lattice artefacts. The 
results for meson masses and the string tension 
have been presented already in section 2. Also 
Bielefeld started to work with Wilson fermions. 
Preliminary results on the phase diagram for tree 
level-improved glue and the clover action were 
presented at this conference Qualitatively, 
the behavior of the kappac and kappat lines is 
very similar to the standard one. 

3.4. Wilson fermions Np ^ 2 

The phase diagram has also been studied for 
the number of flavors differing from 2. All anal- 
yses have been carried out with standard Wilson 
fermions on lattices with Nr = 4 so far. 

The case Np = 1 has been looked at in . For 
heavy quarks the deconfinement transition shows 
first order behavior as a finite size analysis of the 
Polyakov loop reveals. The first-order transition 
weakens when the quark mass is decreased. The 
end-point of the transition is estimated to occur 
at about 1.1 GeV. 



The cases Np = 3 and 4 have been studied in 
Tsukuba |Q . In both cases the phase diagram is 
very similar to the one with 2 flavors, in particu- 
lar the Kc line forms a cusp. Regarding the nature 
of the chiral transition, for Np > 3 one expects 
flrst order in the continuum limit ||20| . At large 
quark masses, away from the cusp, one indeed ob- 
serves first order behavior. When the quark mass 
is lowered however, the Np — 4 data shows a 
weakening of first-order signals. For Np = 3, one 
has to use an approximate algorithm e.g. the Hy- 
brid R. At a time step size of St = 0.01 two-state 
signals were observed also close to Kc- These sig- 
nals weaken when the time discretization errors of 
the algorithm are decreased {St = 0.005). Thus, 
for both Np = 3 and 4 the order of the transition 
is still unclear. 

4. EQUATION OF STATE 

Thermodynamic quantities like pressure or en- 
ergy density receive substantial contributions 
from high momentum modes, p — ttT. Since 
on the lattice these are distorted by UV cut-off 
effects e.g. the energy density in the infinite tem- 
perature limit deviates considerably from the con- 
tinuum Stefan-Boltzmann value. For the pure 
gauge theory in the standard Wilson discretiza- 
tion the corrections are 



^SB 



l + ^f^ 
21 \Nr 



+ - 



0{ 



Nr 



(10) 



This effect is even larger for the fermionic part of 
the energy with standard staggered quarks. 



441 \Nr \Nr ' 



(11) 



On the other hand, the signal for this quantity 
vanishes proportional to ^/N^. Bulk thermody- 
namic quantities were investigated in the stan- 
dard discretization, both quenched j^] and with 
2 fiavors of staggered quarks . While for the 
quenched theory it was possible to extrapolate to 
the continuum, in view of the analytic results. 
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Figure 11. The gluonic part of the energy density 
at high temperature computed on lattices with 
finite temporal extent in units of the continuum 
Stefan-Boltzmann value for various gauge actions 
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Figure 12. The fermionic part of the energy den- 
sity at high temperature computed on lattices 
with finite temporal extent in units of the con- 
tinuum Stefan-Boltzmann value for various stag- 
gered fcrmion actions Pq]. 



Figures ^ 12, this seems to be very difficult with 
dynamical fermions. 

In the infinite temperature, Stefan-Boltzmann 
limit, cut-off effects are reduced at tree level. Per- 
turbative or tadpole improvements do not have 
any effect in this limit. The size of the remaining 
corrections to the Stefan-Boltzmann limit can be 
read off Figures [l^jlj where a variety of different 
actions are compared with each other |4j,^|] . 

The effect of various improvement schemes at 
temperatures as low as the critical temperature 
has been tested in quenched simulations Q . For 
the pressure, it seems that tree level improve- 
ment is the leading effect, although one would 
have expected that close to the transition infra- 
red modes and their improvement would be more 
important. Tadpole improvement seems to have 
an effect though for the interface tension \ ^\ . 

A first attempt to analyze bulk thermodynamic 
quantities by means of an improved fermion dis- 
cretization scheme, the Naik action, has been 
published recently ||l^. At the moment, there 
are investigations underway which try to estimate 
the effect of various improvement strategies on 
the restoration of flavor symmetry [^8|- ^^ . It 
remains to be seen how much this would help to 
extract energy density or pressure closer to the 
continuum limit at finite temperatures. 



5. SCREENING 
MASSES 



LENGTHS 



AND 



The study of hadronic excitations at finite tem- 
perature is an important subject as it reveals 
more details about the properties of the plasma 
phase as well as the approach to the transition 
where temperature dependencies could mimick 
plasma signatures. 

The temperature dependence of pseudoscalar 
and vector meson (screening) masses has been 
studied in |51|. The masses were determined 
at temperatures slightly below and above the 
transition, T ~ 0.9Tc and T ~ 1.2Tc respec- 
tively, and compared with zero temperature re- 
sults. The investigation was carried out with the 
Sheikholeslami-Wohlert fermion action applied to 
quenched gauge configurations. The results, ob- 
tained in the quark mass region ranging from light 
to approximately the charm quark masses do not 
exhibit any strong temperature effect for 0.9Tc. 
At 1.2Tc the data are approaching the high tem- 
perature, free quark limit, M^^p<2ttT. This con- 
firms earlier results for staggered fermions . 

Heavy quarks at temperatures around the criti- 
cal one are investigated in . The heavy quarks 
in the range tUc < m < mt are simulated by 
means of an NRQCD action applied to quenched 
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configurations which have been generated with a 
(1 X 2) tree level improved gauge action. The 
propagation of quarkonia states is followed in the 
time direction. For that purpose, the investiga- 
tion is carried out on anisotropic lattices with a 
large anisotropy ratio ^ = a^jar = 4.65 in order 
to have enough Matsubara frequencies. So far, 
the analysis was done for the ^5*1 ground and first 
excited state. Below Tc, at about O.STc no tem- 
perature effect was seen. At 1.2Tc, the propagator 
which is dominated by the ground state at large 
time separations t from the source becomes flat- 
ter than the zero temperature propagator at large 
t. This effect can be interpreted as a decrease in 
the mass by a small amount of about 12 MeV at 
the lightest quark mass simulated. The effect be- 
comes weaker with increasing quark mass. The 
first excited state, projected on by the same trial 
wave function as at T = 0, undergoes a larger 
change of about - 240 MeV at the charm quark 
mass. At least qualitatively, the results are in ac- 
cord with the expectations from a Debye-screened 
potential model, namely that smaller states feel 
the screening less. 

Investigations of the Debye mass itself have 
been presented by jsj,^. The aim of both 
groups is to confront the lattice data with per- 
turbative results. It is known since long 
that high temperature QCD is plagued by in- 
frared divergencies. These are cured at 0{gT) 
by the Debye or electric gluon mass mg. It is 
widely believed that the magnetic mass m™ is 
non- vanishing and acquires a value of 0{g^T). 
Such a mass would solve the infrared problem 
at this order but cannot be calculated in pertur- 
bation theory. If rrim is non-vanishing next-to- 
leading order perturbation theory predicts pTj for 
the electric mass 



(1 



N 



9- 



6 



-[In 



2m„ 



(12) 



where 



2t: \l 2N + Nf 
rrieQ denotes the leading term TOeq 
JN/3 + Np/GgT. The methods used by [|| and 
pSil are quite different. In [s^ dimensional reduc- 
tion is applied and the 3d effective theory is sim- 
ulated. The Debye mass is then extracted from 
the correlation of a gauge-invariant operator p8| . 



The authors of ||55|] generate SU{2) gauge config- 
urations with the standard as well as an improved 
action. They analyze gauge- variant gluon and 
Polyakov loop correlations after fixing to Landau 
gauge. Their data show a non- vanishing mag- 
netic gluon mass which varies with temperature 
according to mm{T) = cg'^{T)T. For the electric 
mass, the results of the two groups differ some- 
what quantitatively. However, both groups show 
clearly that (next-to-leading order) perturbation 
theory describes the data only at very large tem- 
peratures ofOiW^Tc). 

6. QCD AT FINITE DENSITY 

It is well known that at finite chemical potential 
/i 7^ the effective action of QCD becomes com- 
plex. This prohibits the use of standard Monte 
Carlo methods to evaluate the path integral. A 
way to circumvent these difhculties is provided 
by what is known as the Glasgow method [59[| . 
This method was applied to QCD at intermedi- 
ate and recently at strong coupling [60|. In the 
latter case mean- field computations |61] and re- 
sults from a monomer-dimer simulation | ]62| are 
available. The results of show a finite density 
transition at a critical value /ic which agrees with 
the mean-field as well as with the monomer-dimer 
number. However, also an onset fj,o at which the 
number density starts to deviate strongly from 
is observed. The value of /Xo coincides with 
TO7r/2 for smaller quark masses so that in the chi- 
ral limit the chiral condensate would immediately 
drop from its zero chemical potential value when 
fi is increased from 0. Similar findings were ob- 
tained at intermediate coupling. The same be- 
havior was also seen in the quenched approxima- 
tion where it was traced back to the failure of 
this approximation at finite ^ |l63| . The reason 
for this behavior in full QCD is not clear. One 
might speculate |^ that sampling configurations 
at /i = might not catch the right physics at fi- 
nite /i. Similar observations on /Xo were reported 
in where the modulus of the determinant was 
used in the partition function. A different ap- 
proach has been presented in | ]65[ | where pertur- 
bativcly irrelevant four-fermion interactions were 
added to the QCD action. This way one is able to 
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simulate at vanishing bare quark mass, but it has 
to be clarified yet whether this approach really 
produces only physical results. 

7. CONCLUSIONS 

The results of last year's research on finite 
temperature lattice QCD with dynamical quarks 
have been reviewed. At currently accessible 
quark masses the critical temperature for 2 fla- 
vors is being estimated as around 150 MeV. This 
number is, not unexpectedly, quite smaller than 
the quenched value of 270(5) MeV. Correspond- 
ingly, at fixed Nr the computations with quarks 
have been carried out at considerably larger lat- 
tice spacings compared to quenched simulations. 
Moreover, going to lighter quarks drives one to 
stronger coupling. One may speculate that lat- 
tice spacing effects are the reason for the very 
unexpected behavior at the chiral/deconfinement 
transition for both standard staggered and Wil- 
son fermions while improving only the gauge part 
already seems to lead to remarkable agreement 
with the theoretical predictions. In addition, 
the computation of the energy density of a free 
fcrmion gas shows that the discretization effects 
can be large for fermions when the standard lat- 
tice transcriptions are used. Thus it is highly de- 
sirable to find improved actions for the analysis 
of the thermodynamics with quarks. 
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